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ABSTRACT 

Gravitational lensing of a gamma ray burst (GRB) by a single point mass will 
produce a second, delayed signal. Several authors have discussed using microlensed 
GRBs to probe a possible cosmological population of compact objects. We analyse 
a closely related phenomenon; the effect of microlensing by low to medium optical 
depth in compact objects on the averaged observed light-curve of a sample of GRBs. 
We discuss the cumulative measured flux as a function of time resulting from delays 
due to microlensing by cosmological compact objects. The time-scale and curvature of 
this function describe unique values for the compact object mass and optical depth. 
For GRBs with durations larger than the detector resolution, limits could be placed 
on the mass and optical depth of cosmological compact objects. The method does not 
rely on the separation of lensed bursts from those which are spatially coincident. 
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1 INTRODUCTION 



Press & Gunn (1973) were the first to propose that a cosmo- 
logical abundance of dark compact objects could be detected 
by gravitational lensing of more distant sources. If a point 
mass lies along the observer-source line of sight, the rela- 
tive motion between the lens, source and observer produces 
a change in the magnifications of two lensed images. The 
presence of foreground compact objects is therefore detected 
through a change in the observed flux of the background 
source (termed microlensing) . This effect has been used suc- 
cessfully in the search for compact objects in the halo of the 
Milky Way galaxy (e.g. Alcock et al. 2000). Also, microlens- 
ing due to stars in a galaxy at moderate redshift has been 
observed in the gravitationally lensed quasar Q2237-f 0305 
(Irwin et al. 1989; Corrigan et al. 1991). The short dura- 
tions of gamma ray bursts (GRBs) offer an alternative way 
to study microlensing and hence search for a cosmological 
population of compact objects, through observation of re- 
peating bursts. 

In this paper we assume GRBs to be at cosmological 
distances (Paczynski 1995), and therefore probable sources 
for gravitational lensing (Turner, Ostriker & Gott 1984). 
Paczynski (1986, 1987) noted that while multiple images 
of a single GRB cannot be angularly resolved by present 
day detectors, their relative delay may be longer than the 
burst duration so that the lensed images could be resolved 



temporally (a pair of lensed GRB images that are produced 
by a single mass {Mco) have a relative time delay that is 
At ~ 50sec X [Mco/IQ^Mq) (Mao 1992)). Microlensing of 
a GRB by a compact object is therefore observed as a GRB 
that repeats. The utility of GRBs to explore cosmological 
dark matter, as well as the possible inference of properties 
of the GRB population itself was discussed in detail by Blaes 
& Webster (1992). The short event duration, as well as the 
transparency of the universe to gamma rays make GRBs 
ideal probes of dark matter in the form of compact objects 
over a wide range of masses. 

Microlensing of existing and potential catalogues of 
GRBs have been used to discuss the cosmological abundance 
of compact objects. Marani et al. (1999) use non-detections 
of lensed images from the BATSE and Ulysses catalogues 
to set conservative limits on dark compact objects with 
masses between 10"" and 10"'^ Mq. Also a universe pro- 
posed by Gnedin & Ostriker (1992) with Mco ~ IO'^'^Mq 
and Q.M ~ ^co ~ 0.15 was ruled out at a confidence level 
of 90%. This scenario had been previously investigated in 
detail, using a 3-D lensing code by Mao (1993). He found 
results that did not depart significantly from those obtained 
by a single-screen approximation. Mao (1992) and Gross- 
man & Nowak (1994) estimate that the waiting time for one 
lensed pair due to an intervening galaxy to be observed in 
the BATSE catalogue is between one and 10 years, and con- 
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elude that it is not certain that such a lensed pair will be 
found by BATSE. 

As in quasar lensing, the detection of a lensed pair 
is separated from spatially coincident events by the com- 
parison of light-curves and spectra. However the presence 
of noise and the faintness of images will alter the light 
curves so that macro images appear dissimilar (Wambsganss 
1993; Nowak & Grossman 1994). In addition, if a foreground 
galaxy is responsible for the lensed pair then Williams & Wi- 
jers (1997) find that microlensing by individual stars in that 
galaxy can smear out image light curves, further increas- 
ing the chance of their being mis-classified as two spatially 
coincident events. Light-curve similarity of a lensed pair re- 
quires the source to be much smaller than the Einstein ring 
radius of the lens (projected into the source plane), a con- 
dition likely to be filled in the case of GRBs due to the 
short event duration for all but the smallest lenses. However 
an anisotropic source (eg. a beamed source) is effectively 
viewed at two different angles in the two images. This raises 
the possibility that the time variation of the two images 
may be different (eg. Babul, Paczynski & Spergel 1987), al- 
though Blaes & Webster (1992) find that a source should 
be isotropic to lensing provided that the beaming is not too 
strong. Paczynski (1987) noted that the potential difficul- 
ties in identifying lensed pairs will make confirmation of the 
cosmological origin of the GRBs difficult through a lensing 
argument. On the other hand, if we assume that GRBs have 
a cosmological origin, the mis-classification of events will 
lead to an underestimate of a compact object population. 

The microlensing effect of a cluster of objects on a GRB 
was first discussed by Paczynski (1987). He found that a sin- 
gle instantaneous burst lensed by a screen of objects such as 
that produced by a galactic halo or a cluster of galax;ies re- 
sults in the observation of many repeats of the original burst 
separated by different delays. Examples of the light-curves 
produced can be found in Paczynski (1987) and Williams & 
Wijers (1997). At moderate optical depths the first images 
have magnifications that are not correlated with their delay, 
however the faint images arriving later have observed fiuxes 
that decrease monotonically with time. Williams & Wijers 
(1997) investigated the fiux weighted rms of the delay of 
clusters of microimages and found this to be a sensitive func- 
tion of optical depth and shear. The temporal spread of the 
cluster of microimages was shown to be related to the longest 
side of the area of the image plane containing microimages 
that account for most of the macroimage flux. This scale 
length is considerably larger than that of the separation of 
a lensed pair produced by a single microlens, and so the 
spread in the flux weighted arrival time is also considerably 
larger than the relative delay between images produced by 
a single mass. This phenomenon allows smaller masses to 
be probed modulo a detector resolution and intrinsic GRB 
duration. 

A non-negligible fraction of the baryon content of the 
universe may exist in stellar remnant form (eg. Kerins & 
Carr 1994). In addition, Galactic microlensing searches (e.g. 
Alcock et al. 2000) provide evidence for a signiflcant frac- 
tion of the Galactic halo mass being in the form of stellar 
mass compact objects, an interpretation which is supported 
by observations of high proper motion, cold white dwarfs 
(Ibata et al. 2000; Hodgkin et al. 2000). While intervening 
galaxies produce relatively large delays, GRBs with micro- 



second durations or variability would resolve microlensing 
by cosmological stellar mass objects (Nemiroff, Norris, Bon- 
neU & Marani 1998). Nemiroff et al. (1998) note that there 
is no known fundamental reason for the non-existence of 
micro-second GRBs, and find that a significant rate of these 
might exist and go undetected by current telescopes. They 
find that at a given fiux level there may be more GRBs with 
durations between 1 and 2 milliseconds than between 8.192 
and 16.384 seconds (BATSE duration bin), and an abun- 
dance of microsecond spikes that is an order of magnitude 
greater. However, Nemiroff et al. (1998) also note that the 
converse possibility of short duration spikes not existing at 
all is equally consistent with the analysis of current data. 

We propose that if short duration bursts were to ex- 
ist then they could be used to probe the cosmological stel- 
lar mass compact object population without resolving and 
identifying individual micro-images. Also, in the absence of 
short spikes, the light-curves of GRBs can place upper limits 
on typical compact object mass and optical depth. Rather 
than considering the probability of lensing by a single mass 
on individual GRBs, we consider the microlensing effect on 
an ensemble of GRB light curves of a collection of point 
masses at low to moderate optical depth. For a sample of 
microlensed spike GRBs we look at the average value of the 
cumulative fiux as a function of time. 

This paper presents calculations corresponding to the 
microlensing effect of populations of stellar mass objects on 
the observed light-curves of hypothetical, very-short dura- 
tion GRBs. However, the results are also applicable to the 
presently known population of GRBs, with durations of IDs 
of seconds as probes of compact objects having masses of 
~ lO^Af©. The relevant scaling is pointed out where appro- 
priate. 

The remainder of the paper is presented in 3 parts. Sec- 
tion ^ describes the microlensing models and section ^ pro- 
poses a new method for analysing microlensing in GRBs. A 
brief summary is presented in section 0. 



2 MICROLENSING MODELS 

As a crude approximation to the lensing effect of compact 
objects distributed along the line of sight, we assume a 
screen of point masses distributed randomly in a disc. The 
model does not include a continuous matter distribution. We 
use standard notation for gravitational lensing. The Einstein 
radius of a microlens in the image plane is denoted by 
and when projected into the source plane by rjo. The nor- 
malised convergence or optical depth is denoted by k. The 
normalised lens equation for a field of point masses is 

/ \ N.. ■ 

- / 1-7 \ ^ , (f^ -f) 

Here x and y are the normalised image and source posi- 
tions respectively, 7 is the applied shear, and the Xi and 
m? are the normalised positions and masses of the individ- 
ual microlenses. To construct a microlensed light-curve for 
an instantaneous microlensed GRB, Eqn. ^ is solved for the 
micro-image positions at many points along a predefined 
source line through the inversion technique of Lewis et al. 
(1993) and Witt (1993). The time delay and magnification 
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Figure 1. The fraction of flux tiiat is yet to arrive as a function of time. The cases shown are for mean masses of IMq and IOMq for 
optical depths of k =0.025 (thin lines) and k =0.250 (thick lines). 
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Figure 2. The fraction of flux that is yet to arrive as a function of time normalised to Tgg. Left: The curves correspond to the cases 
shown in Fig. i|. Right: The curves correspond to optical depths k =0.025, 0.050, 0.075, 0.100, 0.125, 0.150, 0.175, 0.200, 0.225, 0.250. 
(thicker lines denote larger k) at a mean mass of IOMq . 



are then determined for each micro-image i, given by 
ATi = 



■(1 + Zd) 



and 

_ 1 

^ \detA{x,)\ 
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^mnn{\x,^x^\)\ (2) 



i=o 



(3) 



where detA(xi) 



dyi I 



dx2 



Here Ds , Dd and D^s are the angular diameter distance be- 
tween the observer and the source, the observer and the 
lens, and the lens and source, c is the speed of light and Zd 
the lens redshift. We note that the delay is proportional to 
{1 + Zd), but is not explicitly dependent on the source red- 
shift. Also, due to the dependence on AT is proportional 
to the mean mass. 
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At a sufficiently large angle from the point source, there 
is a low magnification image located very close to each point 
mass. All solutions of the lens equation must therefore be 
found in a region that contains a sufficient percentage of the 
total macroimage flux. The region of the lens plane in which 
image solutions need to be found to ensure that 99.9% of 
the total macro-image flux is recovered from all points on 
the source line is known as the shooting region. The number 
of stars in the region about any point which collects 99.9% 
of the macro-image flux was calculated by Katz, Balbus & 
Paczynski (1986) and is given by: 



N, = 300 



{m)2 1(1 — k)2 — 7^ 



(4) 



In the absence of shear (7 = 0) these stars are distributed in 
a disc with a radius Rsd = k,/{N, x (m)). The shooting 
region is defined by the union of these discs centred on the 
points {xi = 2/1/(1 — k),X2 = 2/2/(1 — ft)) corresponding to 
all parts of the source line. The minimum number of stars 
(Nmin) required for the model are contained in a disc having 
a radius _R* which covers the shooting region. We have used 
500 stars in each of our models which is larger than Nmin in 
all cases. We assume that the source is stationary with re- 
spect to the microlenses for the duration of the microlensed 
light-curve. We find the light curves for 10 bursts per mi- 
crolens field distributed along a source track of length lOr/o- 
At optical depths between k = 0.025 and 0.25 (at intervals 
of 0.025) we compute the light curves for 10* bursts. 

The most probable scenario for lensing has the lens ly- 
ing at a distance that is a reasonable fraction of that of the 
source (Turner, Ostriker & Gott 1984). In addition, the dis- 
tribution of optical depth with redshift is reasonably sharply 
peaked. With this in mind we stress that our calculation uses 
a single screen to approximate the lensing effect of all masses 
along the line of sight. As an example we place the popula- 
tion of model GRBs at a redshift of Zs — 1 and a screen of 
compact objects at a redshift of = |. For simplicity we 
assume that all compact objects have a unique mass m. 



3 METHOD OF ANALYSIS 

In this section we consider the analysis of a hypothetical 
sample of spike GRBs that have an intrinsic duration smaller 
than 0.1 milliseconds which we take as the detector binsize. 
From our sample of model microlensed light curves we find 
the cumulative flux (summed over all located microimages) 
that has arrived by the end of each bin. One minus the 
cumulative flux divided by the total flux gives the fraction 
of total flux that is yet too arrive (Fi(AT)) for each burst i 
as a function of time. These curves are averaged over many 
bursts: 



F(Ar) 



F,{AT)/Ni^rsts 



(5) 



Eqn. Q describes the average behaviour of the arrival time 
of lensed flux for a set of microlensing parameters k and m. 
Note the assumption that all lines-of-sight have the same 
average optical depth k. From the work on clusters of GRB 
microimages by Paczynski (1987) and Williams & Wijers 
(1997) we expect F{AT) to have the following character- 
istics. Since the microlensing spread of the GRB is related 
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Figure 3. Contour map showing To. 99 in milliseconds (dark 
lines) and FAT„orm=0.i (light lines) as a function of optical depth 
and mass. 



to the size of the shooting region, which is proportional to 
the Einstein radius of the microlenses and therefore to \/rn, 
AT{m) such that F(Ar) = const is proportional to m. Sec- 
ondly, at larger optical depths there are more microimages 
with large magnification. A given fraction of microlensed 
flux therefore arrives later at higher optical depths and so 
F{AT) is larger for increased k. Note that when in a region 
where (1 — k) — 7 ~ 0, the magniflcation becomes very 
large, and flux from the burst will be observed at a com- 
parable level for a time similar to the geometric delay for a 
trajectory at the edge of the region. 

Fig. hi shows the fraction of flux that is yet to arrive as 
a function of time. Cases are shown for mean masses of IMq 
and IOMq, at optical depths of k =0.025 (thin lines) and 
K =0.250 (thick lines). At each optical depth, the larger mass 
produces longer microimage delays, causing a given fraction 
of flux to arrive later on average. In addition, at a flxed mass 
a given fraction of flux arrives later at larger optical depths 
producing a curve that has a smaller initial drop in the flrst 
bin followed by a decline that is more rapid than that at 
smaller optical depth. However at higher k, F{AT) is larger 
for aU AT. 

We define To.gg as the time AT at which 99% of the total 
microlensed fiux from the GRB has arrived. For binsizes 
that are small with respect to the microlensed spread of the 
event, To.gg scales linearly with microlens mass. To.gg is also 
a function of k and provides a natural scaling unit for time. 
The left hand panel in Fig. ^ re-displays the curves from Fig. 
^ with time normalised by To.gg in each case {F{ATnorm) 
where AT„o7-m = AT/T0.99). At the end of each bin the 
normalised curves are independent of mass. The curvature 
therefore provides a probe of the optical depth. The right 
hand panel of Fig. ^ displays F{ATnorm) for optical depths 
between k =0.025 and 0.250 at a microlens mass of IOM0. 
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We suggest the value of F{ATnorm = 0.1) as an indicator of 
optical depth. 

For values of k, and m in the ranges 0.025 ^ k ^ 0.25 
and 0.1 <, m <, 100, we made 1000 simulations of F{AT) 
(each F{AT) was calculated using 100 model microlensed 
GRB light-curves) . Fig. ^ shows average values of the quan- 
tities To. 99 (dark contours) and F{ATnorm = 0.1) (light con- 
tours) over a range of optical depths and microlens masses. 
In regions combining small values of optical depth with mi- 
crolens masses m < IM©, the contours are nearly parallel, 
however for most of the parameter space, values of To. 99 
and F{ATnorm ~ 0.1) describe a unique set of k and m. For 
samples containing 100 GRBs the variance in values of To. 99 
and F{ATnorm = 0.1) are an order of magnitude lower than 
the corresponding mean. 

We have demonstrated that a sample of GRBs could be 
used to measure the quantities To.gg and F{AT„orm ~ 0.1), 
and that this combination corresponds to measurements of 
K and m. However it will not be known whether the ob- 
served spread in the arrival time of GRB flux is due solely 
to microlensing effects or whether it is intrinsic to the source 
(unless the individual micro-images are resolved). To. 99 and 
F{ATnorm = 0.1) will therefore be upper bounds on the val- 
ues for instantaneous bursts, and so will exclude a region of 
mass - optical depth parameter space rather than measure 
probable values. In addition, we note that Fig. ^ presents 
results covering only 3 orders of magnitude in mass. How- 
ever the binsize, mass and To.gg are all linearly related, al- 
lowing scaling of the mass and To.gg by the binsize divided 
by 0.1 milliseconds. Estimates of upper limits for To.gg and 
F(AT„orm ~ 0.1) which will probe ~ W^AIq compact ob- 
jects are therefore possible using current data with GRB 
durations of ~ 10 seconds. 

The approach described requires that the source be in- 
tegrated for long enough to be sure that the entire event 
flux (or some appropriately large fraction) has been accu- 
mulated. We now describe an alternative approach that does 
not require measurement of the total flux. Fig. ^ demon- 
strates that the shape of T(AT) is unique for each com- 
bination of m and k. For example, a curve produced by a 
sample of GRBs lensed by a low optical depth of high masses 
may have the same value F{AT = To.gg) as a sample lensed 
by smaller masses having a higher optical depth. However 
at all AT < To.gg the former will have a smaller value of 
F{AT). We have produced the average function ^'^^(AT) 
for masses and optical depths covering the parameter space 
shown in Fig ^. At each combination of mass and optical 
depth {rritrue and Ktme), we make 1000 mock observations 
of Fi(AT) each calculated from 100 GRBs. We look for the 
values of m and k such that Fa„(AT) best fits each mock 
observation Fi{AT). Our criteria for best fit is to minimise 
the value 

D = max{\F-Fa^\). (6) 

This procedure provides a likelihood for measuring m and k 
given the true values mtme and utrue- 

Plh{m,K\mtrue,Ktrue) (7) 

The construction of Fobs i^T) for an observed set of GRBs 
then provides estimates of m and k through application of 
Bayes' theorem: 



p(n.,.)^iVyp.(n.,.|„.',.')p.™.(-',.Vm'rf^', (8) 

where Ppriorim, k) is the assumed prior probability for m 
and K, and A'' is a normalising constant. 

To demonstrate the statistical uncertainty we have plot- 
ted in Fig. ^ the mean and variance of retrieved values of 
m and k verses their assumed values {ntrue,Tntrue) at fixed 
mass and optical depth respectively (ie. sections through 
Pih)- As in the previous calculation, the figure demonstrates 
that a sample of 100 GRBs is sufficient to obtain a consistent 
result. The systematic bias (which is smaller than the statis- 
tical uncertainty) in the recovered values is due to the finite 
grid of m and k over which Fav and pm are computed. This 
calculation would obtain upper limits on the true values of 
m and n due to possible intrinsic spread of the GRB. Also, 
the binsize and mass are linearly related allowing scaling of 
the case presented here. 



4 DISCUSSION 

If gamma ray bursts have a cosmological origin then they 
provide a unique opportunity to study the cosmological 
abundance of compact objects through the identification of 
lensed pairs. However there are several mechanisms that may 
lead to the mis-classification of lensed pairs as spatially co- 
incident but independent bursts. In addition, the identifica- 
tion of a lensed pair may be ambiguous if the relative delay 
is smaller than the event duration. Upper limits on the con- 
tribution to the mass density of compact objects determined 
from a lack of lensed pairs may therefore be underestimated. 

We have demonstrated that simultaneous upper lim- 
its on the average optical depth and mass of compact ob- 
jects can be estimated by measuring the average cumula- 
tive flux as a function of time for a collection of bursts and 
without resolving the individual lensed images. By selecting 
sub-samples of bursts, speciflc populations of compact ob- 
jects could be examined. For example, a collection of bursts 
that are spatially coincident with fore ground galaxies (like 
GRB 971214 (Kulkarni et al. 1998; Diercks et al. 1999)) wiU 
probe the compact object population in galactic halos. In 
particular, if a population of short duration (micro-second) 
cosmological gamma ray bursts are identified, limits will be 
placed on the abundance of stellar mass compact objects. 
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